Is dark matter made of mirror matter? Evidence from cosmological data. 
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We present new fast numerical simulations of cosmic microwave background and large scale struc- 
ture in the case in which the cosmological dark matter is made entirely or partly of mirror matter. 
We consider scalar adiabatic primordial perturbations at linear scales in a fiat Universe. The speed 
of the simulations allows us for the first time to use Markov Chain Monte Carlo analyses to constrain 
the mirror parameters. A Universe with pure mirror matter can fit very well the observations, equiv- 
alently to the case of an admixture with cold dark matter. In both cases, the analyses show a clear 
indication of the presence of a consistent amount of mirror dark matter, 0.05 < f2 m i rror h 2 < 0.12. 
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The existence of dark matter in the Universe seems to 
be an unavoidable evidence, as confirmed by all the cur- 
rently available astrophysical observations at scales rang- 
ing from cosmological to galactic. At the same time, its 
nature is still completely unknown, and is limited to its 
qualitative behaviour for the process of structure forma- 
tion. Together with Big Bang nucleosynthesis (BBN) , the 
most powerful cosmological tests for dark matter candi- 
dates are the cosmic microwave background (CMB) and 
large scale structure (LSS) power spectra, with their in- 
creasing precisions due to the huge observational efforts 
of several groups. 

Previous analytical and numerical studies on CMB and 
LSS power spectra [H-@] have only limited the parameter 
space of mirror dark matter, without providing a defini- 
tive cosmological answer to its existence. Here we finally 
address this question. 

As suggested many years ago by Lee and Yang Q, 
to suppose the existence of mirror matter is the sim- 
plest way to restore the parity symmetry of the laws 
of nature. Their idea was later developed by other au- 
thors [8MTT| . and a lot of studies were devoted to it, 
showing its compatibility with all the available exper- 
imental and observational constraints (for reviews, see 
Refs. @,|IM1)- In some cases, as the results of direct 
detection experiments fill or the observations of neu- 
tron stars pjKUI, there are interesting suggestions of the 
existence of mirror matter. 

The original idea was to have a parallel hidden (mirror) 
sector of particles which is an exact duplicate of the ob- 
servable sector. In the modern context of gauge theories 
this implies the existence of an exact parity (mirror) sym- 
metry between two particle sectors, that are described by 
the same Lagrangian and coupling constants, and conse- 
quently have the same microphysics, but where ordinary 



particles have left-handed interactions, mirror particles 
have right-handed interactions Thus, they arc sta- 
ble exactly as their ordinary counterparts. 

The ordinary and mirror particles have the same 
masses and obey to the same physical laws, but the three 
non-gravitational interactions act on ordinary and mir- 
ror sectors completely separately, the only link between 
all of them being the gravity. Since mirror baryons do 
not interact with photons, or interact only very weakly, 
the presence of mirror matter is felt mainly by its grav- 
itational effects, which is exactly the definition of "dark 
matter" . 

Hence mirror matter is a stable self-interacting dark 
matter candidate that emerges if one, instead of (or in 
addition to) assuming a symmetry between bosons and 
fermions (supersymmetry), assumes that nature is parity 
symmetric. 

Besides being a viable and powerful candidate for dark 
matter, the increasing interest on mirror matter is due to 
the fact that it provides one of the few potential explana- 
tions for the recent DAMA annual modulation signal [u| , 
together with the results of other direct detection exper- 
iments (CDMS, CoGeNT, CRESST, XENON) p]|. The 
compatibility of this scenario with BBN constraints has 
already been studied [2(| ■ 

Given its consistency with experiments and observa- 
tions, and the unfruitful attempts to prove the existence 
of the other dark matter candidates, scientific commu- 
nity is facing an emergent question: "is mirror matter 
the dark matter of the Universe (or at least a significant 
part of it)?" One possibility to answer this question is to 
look at the cosmological signatures of mirror particles. 

It is worthwhile to note that the presence of the mirror 
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Astrophysical constraints on self interactions of dark matter 
present in literature are valid only for homogeneous distributions 
of dark matter particles, and are therefore not directly applicable 
to the mirror matter case. 
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sector does not introduce any new parameters in particle 
physics (if we neglect the possible weak non-gravitational 
interactions between visible and hidden sectors) . But the 
fact that microphysics is the same in ordinary and mirror 
sectors does not mean that also macroscopic realizations 
should be the same. The different macrophysics is usu- 
ally parametrized in terms of only two "cosmological" 
free parameters: the ratio x of temperatures of the two 
sectors, in terms of temperatures of the ordinary and 
mirror photons in the cosmic background radiation; the 
relative amount ft of mirror baryons compared to the or- 
dinary ones. 
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where T (T"), n b (Cl' h ), and S {S') are respectively 
the ordinary (mirror) photon temperature, cosmologi- 
cal baryon density (normalized, as usual, to the critical 
density of the Universe), and entropy per comoving vol- 
ume [12j |. 

The present energy density contains relativistic (ra- 
diation) component f2 r , non-relativistic (matter) com- 
ponent £l m and the vacuum energy (cosmological term 
or dark energy) density Da- According to the infla- 
tionary paradigm the Universe should be almost flat, 
fi-tot = fim + fir + fiA ~ 1, which agrees well with the 
results on the CMB anisotropy. Now both radiation and 
matter components contain the mirror componentf0, and 
the matter composition of the Universe is expressed in 
general by 
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where the term f^DM includes the contributions of any 
other possible dark matter particles but mirror baryons. 

At the time of BBN the mirror photons 7', electrons 
e and neutrinos u' e T would give a contribution to the 
energetic degrees of freedom equivalent to an effective 
number of extra neutrino families AN^ ~ 6.14x 4 . Cur- 
rent estimates of A7V„ [2l[ correspond to an upper bound 
x < 0.7, and hence at the nucleosynthesis epoch the tem- 
perature of the mirror sector should be smaller than that 
of the ordinary one, T" < T. 

Due to the temperature difference between the two sec- 
tors, the cosmological key epochs take place at different 
redshifts, and in particular they happen in the mirror 
sector before than in the ordinary one @, [l2[ . The rel- 
evant epochs for the cosmic structure formation are re- 
lated to the matter-radiation equality (MRE) z oq , the 
matter-radiation decouplings (MRD) Zdcc and z' dcc due 
to the plasma recombinations in both sectors, and the 



2 Since mirror parity doubles all the ordinary particles, even if 
they are "dark" (i.e., we are not able to detect them now), what- 
ever the form of dark matter made by some exotic ordinary par- 
ticles, there will exist a mirror counterpart. 



photon-baryon equipartitions Zb 7 and 
occurs at the redshift 
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which is always smaller than the value obtained for an 
ordinary Universe. The MRD takes place in every sec- 
tor only after most electrons and protons recombinc into 
neutral hydrogen and the free electron number density 
diminishes, so that the interaction rate of the photons 



drops below the Hubble expansion rate. 
Tdec U P to small corrections, we obtain 

1 + z dcc - + Z d ec) , 



Since T; 
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so that the MRD in the mirror sector occurs earlier than 
in the ordinary one. It has been shown [H, [l2| that, com- 
paring Eqs. (|3|) and (j4j, for x smaller than a typical value 
x eq the mirror photons would decouple yet during the ra- 
diation dominated period, and the evolution of primor- 
dial perturbations in the linear regime is practically iden- 
tical to the standard cold dark matter (CDM) case. Also 
the photon-baryon equipartition happens in the mirror 
sector earlier than in the ordinary one, according to the 
relation 
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Previous analytical and numerical studies on CMB and 
LSS power spectra [THE have only shown, using a 
qualitative comparison with observations, that: (i) for 
low values of mirror temperatures {x < 0.3) all the dark 
matter can be made of mirror baryons; (ii) for high values 
(x > 0.3) mirror baryons can be present as an admixture 
with CDM. 

Now we are finally able to fit the cosmological param- 
eters and obtain their quantitative estimates. 

TABLE I. Adopted fiat priors for the parameters. 



parameter 


lower limit 


upper limit 




0.01 


0.1 


^cdm h 


0.01 


0.8 


X 


0.05 


0.7 


ft 


0.5 


9.0 


100 e s 


0.1 


10 


T 


0.01 


0.8 


n B 


0.7 


1.3 


ln(10 10 7l s ) 


2.7 
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We have modified the publicly available cosmological 
simulation tools CAMB [23| and CosmoMC [H| in or- 
der to include the effects of mirror matter. Since the 
physics of the mirror particles is the same as our parti- 
cles, we have doubled the equations separately in each 
sector, and considered all the particles when describing 
the gravitational interactions. The recombinations are 
computed separately for each sector. The computational 
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TABLE II. 1-cr constraints on the parameters obtained using different dark matter compositions and cosmological tests. 





parameter 


standard 


mirror 


mirror 


mirror±CDM 


mirror±CDM 








CMB 


CMB±LSS 


CMB 


CMB±LSS 


primary 


fib ft 2 


0.02218 ± 0.00053 


0.02184 ± 0.00069 


0.02188 ± 0.00045 


0.02226 ± 0.00058 


0.02184 ± 0.00053 




o yP 1 


0.1157 ±0.0036 








0.055 ±0.016 


0.0661 ± 0.0062 




n s 


0.961 ±0.013 


0.952 ±0.019 


0.950 ±0.011 


0.969 ±0.015 


0.951 ±0.013 




ln(10 10 A s ) 


3.088 ± 0.032 


3.105 ±0.035 


3.099 ± 0.032 


3.114 ±0.044 


3.099 ± 0.031 




100 6l s 


1.0382 ± 0.0026 


1.0376 ± 0.0036 


1.0363 ± 0.0026 


1.0403 ± 0.0035 


1.0366 ± 0.0028 




T 


0.084 ±0.014 


0.080 ± 0.014 


0.079 ±0.013 


0.086 ±0.014 


0.079 ±0.013 


mirror 


X 


- 


0.29 ±0.10 


0.180 ±0.030 


0.39 ±0.11 


0.170 ±0.034 




ft 




5 75 ± 33 


5 62 ± 1 1 

UiUiJ 1 Will 


3.14 ± 0.71 


2 61 ±015 


derived 




0.295 ± 0.020 


0.345 ± 0.041 


0.343 ±0.016 


0.312 ± 0.042 


0.344 ± 0.034 




n A 


0.705 ± 0.020 


0.655 ± 0.041 


0.657 ±0.016 


0.688 ± 0.042 


0.656 ± 0.034 




Z re 


10.4 ± 1.2 


10.3 ±1.3 


10.3 ± 1.2 


10.8 ± 1.2 


10.2 ± 1.2 




h 


0.685 ±0.016 


0.656 ± 0.032 


0.650 ±0.011 


0.690 ± 0.025 


0.651 ±0.023 




age [Gyr] 


13.84 ±0.12 


13.84 ±0.25 


13.95 ±0.11 


13.56 ±0.25 


13.95 ±0.13 




0"8 


0.827 ±0.022 




0.830 ± 0.028 




0.834 ± 0.034 



times of this modified version of CAMB are consider- 
ably increased, but still fast enough to compute the many 
models needed for a Monte Carlo fit in reasonable times. 
Compared with previous numerical simulations [3- 
, we have used an updated estimate of the primor- 
dial chemical composition of mirror particles present in 
Refs. [TU, [25l - [28j . and a more accurate treatment of the 
recombinations of ordinary and mirror particles using the 
numerical code RECFAST [H]. The new models based 
on CAMB and the more accurate treatment of mirror 
BBN are consistent with the previous ones, but there is 
a strong improvement of the computational time, allow- 
ing us now to constrain the parameters. 

We use a Markov Chain Monte Carlo (MCMC) sam- 
pling of the multi-dimensional likelihood as a function of 
model parameters, based on the computations of CMB 
and LSS power spectra obtained with our modified ver- 
sion of CAMB. 

We sample the following eight-dimensional set of cos- 
mological parameters, adopting flat priors on them with 
broad distributions, as shown in Table [I] the baryon and 
cold dark matter densities Vl^ti 2 and Q. c d m h 2 , the relative 
mirror photon temperature x, the relative mirror baryon 
density /3, the ratio of the sound horizon to the angular 
diameter distance at decoupling 9 S , the reionization op- 
tical depth t, the scalar spectral index n s and the scalar 
fluctuation amplitude A s . The upper limit on x is set by 
the aforementioned BBN limit. In addition, we obtain 
constraints on derived parameters: the matter and dark 
energy densities normalized to the critical density fl m and 
the reionization rcdshift z re , the Hubble parameter 
h, the age of the Universe in Gyr, the density fluctua- 
tion amplitude erg at 8/i _1 Mpc. The runs also include 
weak priors on the Hubble parameter, 0.4 < h < 1.0, 
and on the age of the universe, 10 < age(Gyr) < 20. In 
all the computations we assume scalar adiabatic initial 
conditions in a flat Universe (f^tot = 1), a dark energy 
equation of state with w = — 1, masslcss neutrinos, and 
the number of neutrino families of the standard model 



N eS = 3.046. 

We consider two different chemical compositions of 
dark matter: the case pure mirror and the case mixed 
mirror-CDM. In addition, we perform analyses using two 
different configurations: the CMB only and the CMB 
combined with the LSS. The CMB dataset is provided 
by the WMAP7 team [3(|, which measured the acous- 
tic oscillations of the primordial plasma on degree scales 
with cosmic-variance-limited precision. For the LSS, in- 
stead, we include the power spectrum extracted from the 
SDSS-DR7 luminous red galaxy sample [3l[ limited to 
the length scales larger than k ~ 0.2h Mpc -1 to avoid 
non-linear clustering and scale-dependent galaxy biasing 
effects. For comparison, we run also a MCMC chain for a 
standard ACDM cosmology, with the same assumptions 
and priors, to use as a reference model. The results of 
the runs are shown in Table HH where the estimates of the 
parameters and the 1-er confidence intervals are obtained 
by marginalizing the multi-dimensional likelihoods down 
to one dimension. The corresponding best-fit models arc 
shown in Figs. [1] and [21 

Looking at Table [Til we see that the values of the pri- 
mary cosmological parameters, except obviously for the 
CDM density, do not vary significantly between the stan- 
dard model and both the pure and mixed mirror compo- 
sitions, for both kinds of analysis (CMB and CMB±LSS). 
Going to the derived parameters, there is an increase of 
the matter content of the Universe at the expenses of the 
dark energy. This is partly due to a bigger matter density, 
and partly to the decrease of the Hubble parameter, that 
is anyway compatible with the current estimates. For all 
the models, the non-baryonic matter density is between 
5 and 6 times the baryonic density, in accordance with 
common cosmological analyses. But the most interesting 
result is concentrated in the lines constraining the mir- 
ror parameters. Concerning x, the CMB only analysis 
estimates a pure mirror model with x ~ 0.3, around the 
preferred values able to ex plai n the dark matter direct 
detection experiments (l5l . Il6| , while for mixed mirror 
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FIG. 1. CMB and LSS power spectra for best-fit models with FIG. 2. As in Fig. [T] but for models with baryons, mirror 
baryons and mirror matter obtained using CMB only (dotted matter, and cold dark matter, 
line) or both CMB and LSS (dashed line). For comparison 
we show also the standard model fit (solid line). 



this value is even bigger, around x ~ 0.4. The inclusion 
of the LSS lowers significantly these values at x ~ 0.18 
for both pure and mixed mirror, confirming the higher 
sensitivity of the LSS on x already evidenced in previ- 
ous works [l], [lH . We finally look at the most significant 
parameter, j3, which expresses how much, if any, mirror 
matter is present in the Universe. In this case the results 
obtained using CMB alone or combined with the LSS 
are similar. For pure mirror this value is between 5.5 
and 6, clearly showing that cosmology requires a strong 
presence of mirror matter in order to interpret its ob- 
servables. In case of mirror mixed with CDM, the results 
show similar densities of mirror matter and CDM. This is 
a very interesting result, since it means that, when it can 
freely choose between mirror dark matter or collisionless 
massive WIMPs, cosmology unequivocally requires the 
presence of a consistent fraction of mirror matter. 

In Figs.Q]and[2]thc agreement of the best fits with the 
data is shown, together with the comparison with the 



reference standard model. For the CMB the models ob- 
tained fitting both CMB and LSS data arc almost indis- 
tinguishable, and few differences arc present in the LSS 
power spectra. The exceptions are the models computed 
using CMB only, for which the high values obtained for 
x compromise, as expected, their compatibility with the 
matter power spectrum. 

To summarize, in this work we have obtained two main 
results. First of all, for the first time the two parameters 
describing the mirror matter are constrained. Consid- 
ering the most stringent analyses performed using both 
the CMB and LSS, we obtained x = 0.180 ± 0.030 and 
/3 = 5.62 ± 0.11 for pure mirror, x = 0.170 ± 0.034 and 
j3 = 2.61 ± 0.15 for mixed mirror-CDM. These values 
lie in the range required for interesting consequences on 
observations and experiments. Secondly, but even more 
important, the analyses show a clear indication of the 
presence of consistent amounts of mirror dark matter in 
the Universe. 
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